      subroutine initgeom

      use param
      use allocs

#include <global.h>
      include 'geom.h'
      include 'xsec.h'

! geom.h defaults
	data bcref/0,0,0,0/
	data bcflux0/1,1,1,1/
	data bcin0_p3/0.5,-0.375,-0.125,0.875/
      data bcin0_p1/0.5,0.5/
      data bcm_p1/0,0/
! ls.h defaults
      data idm/1,3,2,4/  !idm(1,1)=1, idm(2,1)=3, idm(1,2)=2, idm(2,2)=4

      end subroutine
